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Abstract 

We present preliminary measurements of branching fractions and CP-violating charge asymmetries 
in i?-meson decays to pvr. The data sample comprises 89 million T(4S') — > BB decays collected 
with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We find the 
charge- averaged branching fractions B{B~^ P^^*^) = (H-O =b 1.9(stat.) it 1.9(syst.)) x 10~^ and 
B{B+ p°Tr+) = (9.3 ± 1.0(stat.) ± 0.8(syst.)) x 10~^; we set a 90% confidence-level upper limit 
of B{B^ — > p^TT^) < 2.5 X 10^^. We measure the CP-violating charge asymmetries A^^p = 
0.23 ± 0.16(stat.) ± 0.06(syst.) and A^^p^ = -0.17 ± 0.11(stat.) ± 0.02(syst.). 
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1 Introduction 



The study of B meson decays into charmless hadronic final states plays an important role in the 
understanding of the phenomenon of CP violation in the B system. Recently, the BABAR experiment 
has performed a search for CP-violating asymmetries in neutral B decays to p^Tr^ [1], where the 
mixing-induced CP asymmetry is related to the angle a = arg [— ^d^ft/^d^l)] Unitarity 
Triangle. However, in contrast to the theoretically clean determination of sin 2/3 in the decay to 
charmonium such as S*^ J/ipK^ [2, 3], the extraction of a from p^ir"^ is complicated by the 
interference of decay amplitudes with different weak phases. Various strategies to overcome this 
problem have been proposed in the literature. One such method is an SU(2) isospin analysis of 
the pTT final states [4]. In the limit of isospin symmetry, the five decay amplitudes B^ P^t^~ , 
B^ P 7r+, B^ P^T^^, p^TT^ and B^ p^ir^ form a pentagon in the complex plane. 

Combining measurements of all the decay rates, mixing-induced and direct CP asymmetries in the 
neutral B modes, as well as charge asymmetries in the charged B modes, allows a determination 
of the phase a that is free of hadronic uncertainties. 

In this letter, we present preliminary measurements of the branching fractions of the decay 
modes^ B~^ — > p^vr"'", B~^ — > p'^tt^, and perform a search for the decay B^ — > p^vr". For the charged 
modes we also present measurements of the CP- violating charge asymmetry Acp, defined by 

r(B- ^ /) - TjB* ^ /-) 

^'^ r{B- ^ /) + r(B+ ^ /) ^ ' 

where T{B^ f) and T{B~^ /) are the B~ and B^ decay rates, respectively. A non-zero Aqp 
requires the presence of at least two amplitudes with different weak and strong phases. 

The measurements of the branching fractions and charge asymmetries use events collected by 
the BABAR detector at the PEP-II asymmetric-energy B Factory. The sample used for the charged 
modes consists of 88.5 x 10^ BB pairs collected at the T{4S) resonance ( "on-resonance" ) , while 
the B^ —>■ p^tt" search uses a slightly larger sample of 88.9 x 10^ BB pairs. We use an integrated 
luminosity of 9.6 fb~^ collected approximately 40MeV below the T{4S) ("off-resonance") for back- 
ground studies. We determine the yields and charge asymmetries using a maximum likelihood 
(ML) fit. 



2 The BABAR Detector 

A detailed description of the BABAR detector can be found in Ref. [5]. Charged-particle momenta 
are measured in a tracking system consisting of a 5-layer double-sided silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) filled with a gas mixture based on helium and isobutane. The 
SVT and DCH operate within a 1.5 T superconducting solenoidal magnet. The typical decay vertex 
resolution is around 65 pm along the beam direction for the fully-reconstructed Bp^ (referred to as 
signal B hereafter), and around 185 pm for the inclusively-reconstructed rest of the event. Photons 
are detected in an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals arranged 
in barrel and forward end-cap sub-detectors. The average tt^ mass resolution is 7MeV/c^. The 
flux return for the solenoid is composed of multiple layers of iron and resistive plate chambers for 
the identification of muons and long-lived neutral hadrons. Tracks are identified as pions, kaons 
or protons by the Cherenkov angle 9c measured with a detector of internally reflected Cherenkov 

'^If not stated otherwise, charge conjugation is implied throughout this document. 
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light (DIRC). The typical separation between pions and kaons varies from 8a at 2GeV/c to 2.5a" 
at 4GeV/c, where a is the average 9c resolution. Lower momentum kaons are identified with a 
combination of 9c (for momenta down to 0.7GeV/c) and measurements of energy loss, dE/dx, in 
the DCH and SVT. 

3 Discriminating Variables 

To reject background from continuum e^e~ qq {q = u, d, s, c) events and other B decays, we use 
the following discriminating variables: 

• ruES- the beam-energy-substituted mass is defined by 



where s is the square of the center-of-mass (CM) energy, Ei and pi are the total energy 
and three-momentum of the initial e~^e~ state in the laboratory frame, and pb is the three- 
momentum of the B candidate in the same frame. Signal events populate the niES region 
near the B mass. Background events usually have a wider or even flat uies distribution. 

• AE: the difference between the reconstructed energy of the signal B candidate and the 
beam energy in the CM frame. The AE distribution for signal events peaks around zero. 
Backgrounds from other i?-meson decays peak at different AE' depending on the number of 
charged and neutral particles in the decay. Two-body and four-body decays populate the 
positive and negative regions of AE, respectively, while three-body decays have AE values 
near zero. The AE is shifted if a misidentified charged track is used to reconstruct the signal 



• Output of a multivariate analyzer : the dominant background for charmless B decays are 
continuum events. To enhance discrimination between signal and contimium, three neu- 
ral networks (NN), one for each signal mode, arc constructed. Each NN is trained with 
off-resonance data and simulated events of its corresponding signal. All NNs use: the recon- 
structed p mass; cos^ttj where 9,^ is the angle between the momentum of the 7r"'"(7r°) from the 
p^{p^) decay and the B momentum in the p rest frame; the event shape variables Lq = ^iP* 
and L2 = J2iPi ^ ! '^o^(^Tb i)!^' where p,* is the CM momentum of the track or neutral object 
i, belonging to the rest of the event, and 9^^ ^ is the angle between the momentum of track i 
and the B thrust axis Tb in the CM frame; | cos{9b,z)\, the cosine of the angle between the 
B momentum and the z axis (along the beam direction) in the CM frame, and | cos(^Tb,2)|) 
the cosine of the angle between the B thrust axis and the z axis in the CM frame. 

4 Event Selection and Reconstruction 

Signal B candidates are reconstructed from combinations of three-pion final states that must be 
7r+7r°7r° for B+ p+n^, 7r+7r-7r+ for B+ p°7r+, and Tr+TT~TT° for B° p^ir^. The charged 
tracks must be inconsistent with being an electron based on dE/dx measurements, shower shape 
criteria in the EMC, and the ratio of shower energy and track momentum. Charged tracks must 
also be inconsistent with being a kaon or a proton. The photon candidates used to form a 7r° must 
have an energy greater than 50 MeV in the laboratory frame, and exhibit a lateral shower profile 




(2) 



B. 
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of energy deposition in the EMC consistent with an electromagnetic shower. The invariant mass 
777.(77) of the photon pair must satisfy 0.11 < 777,(77) < 0.16GeV/c^. 

Two of the three final state pions are used to form a charged or a neutral p candidate (the 
third pion is referred to as the bachelor pion hereafter). For the p'^ir^ mode, the p'^ 

candidate is reconstructed from the positively charged track and one of the two tt" candidates. For 
the p^tt"*" and — > p°7r° modes, the candidate is reconstructed from two oppositely- 

charged tracks. The mass of the p candidates must satisfy 0.4 < m{TT'^'K^) < 1.3GeV/c^ for 
B+ p+7r°, 0.53 < m(7r+7r-) < O.QGeV/c^ for B+ p°7r+, and 0.4 < ^(vr+vr-) < O.gCcV/c^ for 
B^ p^ir^. The B^ p^n^ mass cut has been tightened to remove S-related backgrounds such 
as B+ ^ fo{980)TT+{K+), K^n+. The ^ /9°7r° mass cut is tight to remove S° ^ /o(980)7r°. 
For the B^ — > p^tt^ mode, the corners of the Dalitz plot where p^tt^ interferes with the dominant 
p^TT^ are removed by requiring that the invariant masses 77i(7r+7r^) and 7n(7r~7r'') be greater than 
1.3 GeV/ . To take advantage of the helicity structure oi B pir decays, we require | cos 6t^\ > 0.25. 

To reject background from p^K'^ and B^ p~K^ decays, only p^ir'^ candidates with 

bachelor tracks within the geometrical acceptance of the DIRC are considered. The number of 
photons in the DIRC that are associated with the bachelor track must not be less than 5. 

For the B~^ p'^ir^ and B^ — > p^Tr" modes, the invariant mass of either track and the vr'^ must 
be less than 5.14 GeV/ to reject two-body B background. 

B candidates are required to satisfy kinematic fit-region cuts. For S+ p^ir^ decays, can- 
didates must satisfy 5.20 < ruES < 5.29 GeV/ and -0.15 < A£; < 0.10 GeV. B+ p^Tx+ 
candidates must satisfy 5.23 < uies < 5.29 GeV/ and -0.05 < < 0.05 GeV. The B+ p°7r+ 
analysis benefits from a better A£' resolution due to the absence of neutral pions; the tight cut on 
A£' helps to remove four-body B background more effectively. B^ p^n^ candidates must satisfy 
5.23 < m,ES < 5.29 GeV/ and -0.15 < A£; < 0.10 GeV. 

For the B^ — ^ p^TT^ mode, we remove background from charmed decays B ^ D X, D — >■ 
K~^n~ or 7r"'"7r", by requiring that all pairs of oppositely-charged tracks have invariant masses either 
smaller than 1.844 GeV/ or greater than 1.884 GeV/ c^, assuming both kaon and pion hypotheses 
for the positively-charged track. 

The final samples of signal candidates are selected with a cut on the NN outputs for all three 
decay modes. For example, the NN cut for the B~^ — > p^Tr~^ decay mode retains 85% (11%) of the 
signal (continuum) events. 

In each event, final state particles other than the three pions that form the signal B meson are 
assumed to belong to the other B meson. These particles are used to tag the flavor of the other 
B meson and to inclusively-reconstruct its vertex for decay time determination. In this letter, this 
other B is referred to as -Stag- 

For the B^ p^-K^ mode, we use the proper decay time as a discriminating variable in the 
ML fit. The time difference At is obtained from the measured distance between the z positions 
(along the beam direction) of the BpO^^o and Stag decay vertices, and the known boost of the e+e~ 
system. The vertex of the Stag is reconstructed from all tracks in the event except those from 
the BpO^o , and an iterative procedure [2] is used to remove tracks with a large contribution to the 
vertex x^. An additional constraint is obtained from the three-momentum and vertex position of 
the BpO^o candidate, and the average e+e^ interaction point and boost. We require \At\ < 20 ps 
and a{At) < 2.5 ps, where a{At) is the error on At estimated on an event-by-event basis. 

Approximately 33% (7% and 8%) of the events from signal Monte Carlo (MC) simulation have 
more than one candidate satisfying the selection in the B~^ p~^TT^ {B~^ p^-K^ and B^ — > p^ir^') 
decay mode. In this case, we choose the candidate with the reconstructed p invariant mass closest 
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Table 1: Signal efficiencies (e), fractions of misreconstructed signal events (fscf), and fractions 
of misreconstructed signal events with wrong B candidate charge (ujq) in selected MC-simulated 
events. The last row gives the numbers of selected on-resonance events entering the maximum 
likelihood fits. 

e [%] 17.6 ±0.1 29.0 ±0.1 20.2 ±0.1 

/scf [%] 38.6 ±0.2 7.1 ±0.1 9.4 ±0.2 

coQ [%] 8.1 ±0.1 1.6 ±0.1 

Selected events 13104 8498 6648 



to the nominal p mass [6] . Any chosen candidate from a signal event that is reconstructed from one 
or more wrong particles (charged or neutral) is hereafter referred to as misreconstructed signal. 

The signal selection efficiencies are determined by applying the selection criteria to MC. Effi- 
ciencies (including misreconstructed signal) of the three decay modes are summarized in Table 1. 
Also given in Table 1 are the fractions of misreconstructed events in the selected signal samples. 
Misreconstruction is mostly due to p candidates that include a random low-momentum pion. Some 
of the misreconstructed events are assigned an incorrect B charge in charged B decays. The total 
number of events that enter the likelihood fits are given in the last row of Table 1. 

While flavor information of the i?tag is not explicitly used in our analyses, event properties used 
to categorize flavor-tagged events can be exploited to distinguish signal and (mainly) continuum 
background because they populate among tagging categories in dramatically different ways. Five 
mutually exclusive tagging categories are defined: Lepton, Kaon, NTl, NT2, and Untagged. They 
are determined by a tagging algorithm [2] relying on the correlation between the flavor of the 
b quark and the charge of the remaining tracks in the event, after removal of the tracks from the 
B ^ pn. The fractions of signal events in each tagging category are measured from data using a 
control sample of fully-reconstructed B decays [2] . Tagging fractions for B backgrounds are taken 
from Monte Carlo. Separate continuum background yields for each category are free to vary in the 
maximum likelihood fit. Using tagging categories, we decrease the statistical error in the signal 
yield by about 10%. 

5 S-related Backgrounds 

We use MC simulation to study the cross-feed from two-body, three-body and four-body charmless 
B decays, as well as from inclusive h ^ c decays. The branching fractions of unmeasured decay 
channels are estimated within conservative error ranges. The charmless modes are grouped into 
eighteen (i?+ —> p+TT^), fifteen (i?"*" — > p^Tr"*") and seventeen {B^ p^ir^) classes with similar kine- 
matic and topological properties. Two additional classes account for the neutral and charged 6 — ^ c 
decays. For each of the background classes, a component is introduced into the ML fit, with a fixed 
number of events. Contributions to the systematic error by each B background mode is obtained 
by varying its yield by one standard deviation. For unknown modes with only estimated branching 
ratios, the uncertainties are divided by s/S and taken as the standard deviations. Tables 2, 3 and 
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4 summarize the dominant B background modes to P'^t^^, —>■ p^vr"'" and B^ —>■ p^ir^, 

respectively. 

6 The Likelihood 

We use unbinned extended maximum likelihood fits to extract the pir event yields, the charge 
asymmetries, and other parameters used to model signal and background events. The fits minimize 
the quantity — 21n>C, where £. is the total likelihood defined over all tagging categories k by 

£ = ne-^^nA,fe, (3) 

k i=l 

and where ATj.. is the sum of the signal and continuum yields (to be determined by the fit) and the 
fixed 5-background yields, are the numbers of observed events in category k, and Ci^k is the 
likelihood computed for event i. Note that no tagging information is used in the /9°7r+ fit. The 

data sample of each mode is assumed to consist of signal, continuum background and B-background 
components. The variables uiesi and the NN output discriminate signal from background. 
For B^ — > p^TT^ ^ the variable Ai is used to obtain additional background discrimination. 

The likelihood Ci^k for event i is the sum of the probability density functions (PDF) of all 
components, weighted by the expected yields for each component, 

Nb 

C^, k = N^e^V^^ + NfVl% + E , , (4) 

where 

• is the number of signal events in the entire sample. For the p'^ir^ and p^ir'^ modes, the 
charge asymmetries are introduced by multiplying the signal yields by ^(1 — QbAcp), where 

is the charge of B candidate. 

• Ck is the fraction of signal events that are tagged in category k. 

• N^'^ is the expected number of continuum background events that are tagged in category k. 

• The PDFs , V^^' and the likelihood terms C^j^ arc the product of the PDFs of the discrim- 
inating variables. The signal PDFs are thus given 

P^^^^+(NN) for the charged B decay modes, and by vf''° = VP°''\mEs) • Vf°''\AE) ■ 
rf^'^NN) ■ v('"°{At) for B^ p'^ir^. The PDF for continuum events is denoted Vf^^. 
The likelihood term Cf^ j. corresponds to the jth BB background's contribution of the Nb 
B background categories. Correlations between the variables are usually neglected except 
that for B backgrounds we use two-dimensional PDFs for m es and AE to model the sizable 
correlations. 

The signal PDFs are decomposed into two parts with distinct distributions: signal events that are 
correctly reconstructed and misreconstructed signal events. Moreover, for the charged B modes 
we distinguish misreconstructed signal events with right-sign B charge from those with wrong-sign 
B charge in the likelihood. Their individual fractions are taken from MC simulation and given in 
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Tabic 2: B background modes considered in the B~^ — > p'^-K^ maximum likelihood fit. The second 
column gives the branching fractions used ( estimated branching fractions are indicated by an aster- 
isk), and the third column quotes the expected number of events entering into the sample, scaled to 
81.5 fb~^ ("SS.Sx 10^ r{4S) BB decays). The labels "long." and "tran." refer to the longitudinal 
and transverse polarization of the final states, respectively, in B decays into vector-vector mesons. 
K^^*^ refers to the higher kaonic resonances including Kq(1430), iir|(1430) and K*{1680). 
indicates inclusive charmed decays. 



Mode 


BR [10-6] 


N 

-■-'exp 




1.6 ±1.6 


4.8 ±4.8 


B+ ^ ;r+7rO 


5.2 ±0.8 


5.6 ±0.9 


B+ -> X+^o 


12.7 ± 1.2 


4.9 ±0.5 




22.6 ±2.8 


44.6 ± 5.6 


JD — > p I\ 


1 .6 ± i.o 


z.O ± U.D 


B+ -> i^O(7rO^O)^+ 


4.1 ±0.4 


5.4 ±0.5 


^0 ^ A-*"(K+^-)7rO 


8.7 ±5.0 


3.8 ±2.2 




7.5 ± 5.0 


3.9 ±2.6 


B+ ^ i^*+(ii:+7rO)7rO 


4.4 ±2.5 


6.1 ±3.5 


B^ /9+7 


2.3 ± 2.3 


1.3 ± 1.3 


B'^ p'^ p~ long. 




64.8t^^ 


B^ p^p~ tran. 


40.0t|°* 


1.6t?-° 


5+ p'^ p^ long. 




17.0tti 


B+ ^ a+vfO 


35 ± 35* 


25.0 ± 25.0 


B+ ^ {kPtt)+ 


40 ± 26* 


15 ±9.8 


B+ ^ (i^f 


20 ± 20* 


1.2 ± 1.2 




72 ± 54* 


23 ± 17.3 


50 ^ (ii:J*)p)0 


20 ± 20* 


6.0 ±6.0 


Total charmless background 




236.5 ± 79.3 






72.0 ± 16.0 






133.0 ± 30.0 


Total 5 background 




442 ± 86 
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Table 3: B background modes considered in the B~^ p^-K~^ maximum likelihood fi,t. The second 
column gives the branching fractions used ( estimated branching fractions are indicated by an aster- 
isk), and the third column quotes the expected number of events entering into the sample, scaled to 
Sl.hfb'^ (88.5 X 10^ r(45) BB decays). 



Mode 


BR [10-6] 


N 


B^ p-K+ 
B+ ^ p^K+ 
B+ ^ /o(980)if+ 
B+ if0(7r+7r-)7r+ 
B^ ^ K*+{K+ 1x^)71- 
B+ —>■ if*°fif+7r~W+ 

j—^ ' J. \. yj.\. li 1 1\ 

B^ — ^ 7r+a;(7r"^7r~) 
B° — > p'^p~ long. 

B^ p^p^ long. 
B+ r]'{p'^-f)7T+ 
5° ^ a+TT- 
S+ ^ (if 
^0 ^ (if (;*)7r)0 


22.6 ±2.8 
7.3 ± 1.8 
3.9 ± 1.2 

11.7 ±4.0 
9.0 ±0.9 
8.7 ±5.0 
10 3 ± 2 6 
0.14 ± 0.04 

3.5 ± 3.5* 

on 1+8.3 
"^"-'•■'--9.9 

3.0 ± 2.0* 

35 ± 35* 
40 ± 26* 
72 ± 54* 


29.3 ±5.2 
1.1 ±0.3 
4.9 ± 1.5 

1.1 ±0.4 

5.3 ±0.5 

1.4 ±0.8 
11.1 ± 2.8 

3.6 ±1.0 

6 3+^-^ 
"•^-5.5 

1.7 ± 1.7 

7 q+2-2 
' -^-2.6 

2.2 ± 1.5 

5.3 ±5.3 
2.9 ± 1.9 

7.4 ±5.5 


Total charmless background 




91.5±11.3 






19.2 ±5.8 
54.1 ± 13.1 


Total Z? liacliground 




165 ± 19 



Tabic 1. The m^^Si Ai?, and NN output PDFs for signal and B background are taken from the 
simulation except for the means of the signal PDFs for niES and AE in B~^ P^'^'^, which are 
free to vary in the fit. 

The parameterizations of PDFs are discussed in the following. 

• mES- The distribution of correctly reconstructed signal is parametrized using a Crystal Ball 
function [7]. 

The continuum background is described by an ARGUS shape function [8] with floating shape 
parameter ^. 

• AE. We use a Crystal Ball [7] function (sum of two Gaussians) for correctly reconstructed 
signal in p+Tr'^ {p^7r~^ and p^ir^), and the sum of two Gaussians for misreconstructed signal 
events. 
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Table 4: B background modes considered in the B" p^n maximum likelihood fit. The second 
column gives the branching fractions used ( estimated branching fractions are indicated by an aster- 
isk), and the third column quotes the expected number of events entering into the sample, scaled to 
81.9 fb-^ (88.9 X 10^ r{4S) BB decays). 



Mode 


BR [10-6] 


N 


B+ ^ i^+TfO 
B p TT^ 

p-K+ 

B+ ^ p+TfO 

5° ^ ifO(7r+7r-)7rO 

5° ^ /o(980)7r0 

B^ — 7r+7r~7r° (non-res.) 

BO ^ K*^{{KTr)^)j 
B^ p^p~ long. 
— * long. 

S+ ^a^((p7r)+)7rO 
50 ^r/'(p%)7rO 
B+ ^ (i^j;*)7r) + 


5.2 ± 0.8 
12.7 ± 1.2 
ZZ.b ± Z.8 

7.3 ± 1.8 
15.0 ^[l* 
3.5 ± 0.5 
0.0 ± 3.0* 
0.0 ± 5.0* 
8.7 ± 5.0* 
40.2 ± 2.7 

40.0 1:^°* 

30.1 tli* 

35.0 ± 25.0* 
0.0 ± 1.0* 

40 ± 26* 
72 ± 54* 


1.4 ± 0.2 

0.9 ± 0.1 

16. y ± 1.4 

0.7 ± 0.2 

15.8 tin 

5.0 ± 0.7 

+^-^ 
"•^ -0.0 

0-0 

13.6 ± 7.8 

1.7 ± 0.1 

10.9 ™ 
15.8 tti 

5.2 it 3.7 

+2 ° 

-0.0 

2.3 ± 1.5 
3.2 ± 2.4 


Total charmless background 




93.4 it 23.9 






23.8 ± 7.1 
35.2 it 10.6 


Total B background 




152 ± 27 
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Continuum background is modeled by a linear function with a slope that is free to vary in 
the fit. 

• NN output. PDFs for correctly reconstructed and for misreconstructed signal events are 
taken from MC simulation and parameterized using empirical shape-fitting techniques [9]. 
A small discrepancy between data and MC is observed for the NN output distributions of 
control samples using fully-reconstructed — > D~ decays. This is accounted for in the 

systematic error evaluation. 

For B+ p+7r° and B+ p°7r+ /0°7r°), the PDFs describing the NN output for con- 

tinuum events are parametrized by a third-order (fifth-order) polynomial with its parameters 
determined in the fit. 

• At is used in the B p°7r° fit to improve the discrimination against continuum back- 
ground. The distributions for correctly reconstructed signal, misreconstructed signal and B 
background events are treated as decays of neutral or charged S's, convoluted with a At 
resolution function which is the sum of three Gaussians with parameters determined from a 
fit to fully reconstructed B decays [2] . This treatment does not introduce a bias in the signal 
yield according to MC studies. 

The continuum At distribution is parameterized as the sum of three Gaussian distributions 
with common mean, two relative fractions, and three distinct widths that scale the At event- 
by-event error, yielding six free parameters. 

The shapes of the B background PDFs are obtained from MC simulation and parameterized 
using empirical shape-fitting techniques [9]. 

We perform fits on large MC samples with the measured proportions of signals and continuum 
and B backgrounds. Biases observed in these tests are due to imperfections in the likelihood model, 
e.g., unaccounted correlations between the discriminating variables of the signal and B background 
PDFs. The observed signal yields are corrected for these fit biases and the full correction is assigned 
as a systematic uncertainty. 

7 Preliminary Results 

We obtain the event yields 170.8 ± 28.7(stat.) for p+7r°, 232.5 ± 26.4(stat.) for p'^n^ and 15.6 ± 
11.7(stat.) for p^ir^. Assuming equal branching fractions for T(4S') decays into neutral and charged 
B mesons, the yields translate into the branching fractions 

B{B+ -> p+7r°) = (11.0 ± 1.9 (stat.) ± 1.9 (syst.)) x 10*^ , 
B{B+ p°7r+) = (9.3 ± 1.0 (stat.) ± 0.8 (syst.)) x 10^^ , 
B{B^ p°7r°) < 2.5 X 10"^ at 90% C.L. , 

where the first errors are statistical and the second systematic. We find the charge asymmetries: 

AP^p° = 0.23 ± 0.16 (stat.) ± 0.06 (syst.), 
A^^p'^ = -0.17 ± 0.11 (stat.) ± 0.02 (syst.). 

Figure 1 shows distributions of m^^, AE, the NN output and the p mass for S+ — p°7r+, enhanced 
in signal content by cuts on the signal-to-continuum likelihood ratios of the other discriminating 
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Figure 1: Distributions of niES (upper left), AE (upper right), NN output (lower left) and the 
p mass (lower right) for samples enhanced in p^7r~^ signal content using cuts on the signal-to- 
continuum likelihood ratio. The solid curves represent projections of the fit result. The dashed 
curves represent the contribution from continuum events, and the dotted lines indicate the combined 
contributions from continuum events and B backgrounds. For the p mass distribution, the fit has 
been repeated with p-related information removed from the NN. 

variables. The plots of mes, AE and NN correspond to the fit reported here, while the plot of the 
p mass is obtained from a fit with p-related information removed from the NN. 

The statistical significance of the previously unobserved — > p+vrO si gnal amounts to 9.4a, 
which reduces to 6.6a when also considering systematic errors. Figure 2 shows the corresponding 
signal-enhanced distributions of thes and AE. 

For S° /0°7r°, a 90% confidence-level upper limit of 33.2 is obtained on the signal yield using 
a limit setting procedure similar to Rcf. [10]. To obtain an upper limit for the branching ratio, the 
upper limit on the signal yield is shifted upwards by one sigma of the systematic error on the yield, 
and the efficiency and other scaling factors are shifted downwards by one sigma of their systematic 
errors. Figure 3 shows the corresponding signal-enhanced distributions of tues and AE. 

All results are given in Tables 5 and 6 together with the systematic uncertainties discussed 
below. 
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Figure 2: Distributions ofniES (^^fi) o-nd AE (right) for samples enhanced in p'^ir^ signal content 
using cuts on the signal-to- continuum likelihood ratio. The solid curves represent projections of the 

fit result. The dashed curves represent the contribution from continuum events, and the dotted lines 
indicate the combined contributions from continuum events and B backgrounds. 





GeV/c 



AE 



GeV 



Figure 3: Distributions of ruES fl^ft) o-^d AE (right) for samples enhanced in p°7r° signal content 
using cuts on the signal-to- continuum likelihood ratio. The solid curves represent projections of the 
fit result. The dashed curves represent the contribution from continuum events, and the dotted lines 
indicate the combined contributions from continuum events and B backgrounds. 
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8 Systematic Uncertainties 



The systematic errors in the branching fractions arc obtained by adding in quadrature the sys- 
tematic uncertainties in the signal yields, the systematic uncertainties in efficiencies of tracking, 
particle identification, vr'^ reconstruction and the systematic uncertainties on the selection cuts. The 
systematic errors on the Acp measurements are introduced by the uncertainties in the treatment 
of the B background and by possible charge biases of the detector. 

The basis for evaluating the systematic uncertainties on the cuts that are applied in the selec- 
tion process is the differences in AE, niES and NN between on-resonance data and Monte Carlo 
simulation. The differences between data and Monte Carlo distributions of AE and niES are ex- 
tracted from various fully-reconstructed B control samples for the three decay modes. The number 
of DIRC photons cut for the bachelor track in B'^ decay mode will cause 1.0% uncertainty 

on the signal yield. The corrections and uncertainties on the signal efficiencies are summarized in 
Table 5 for the three decay modes. 

We evaluate the systematic uncertainties due to the signal niES-, and NN PDFs with a large 
B data control sample. The small differences observed in the distribution shapes for Monte Carlo 
events and the distribution shapes obtained from the data control sample are used to estimate the 
systematic uncertainty on the signal ruES and AE PDFs. The uncertainties due to the estimated 
fractions of misreconstructed events are obtained from a control sample of fully-reconstructed B 

decays as in Ref. [1]. We perform fits on the large MC samples with the measured proportions 
of p^vr'', p'^7r+ and p^ir^ signals, and continuum and B background. Fit biases observed in MC fits 
are added in quadrature and assigned as a systematic uncertainty of the fit procedure, referred to 
as "fitting procedure" in Table 5. 

The expected yields from the background modes are varied according to the uncertainties in the 
measured or estimated branching fractions indicated in Tables 2, 3 and 4 for the B~^ — > p"'~7r°, B~^ — > 
p'^7r+ and B^ p^ir^ decay modes, respectively. Since B background modes may exhibit direct 
CP violation, the corresponding parameters are varied within their physical ranges. Contributions 
from non-resonant B~^ Tr+vr^vr'^ for the p+vr^ mode and B^ 'k~^it~it~^ for the p^ii^ mode 
are negligible according to our dedicated studies. To check for these types of B backgrounds, a 
fit without p° mass and /9° helicity information in its NN training is performed, and the results 
are compatible with the fit results reported here. The systematic error on the p^vr'^ yield due to 
non-resonant B^ 'k~^it~it^ is considered as part of the B background one, based on Ref. [11]. 

For the B~^ — > p^ir'^ and B^ p^ir^ decay modes, systematic uncertainties due to possible 
interference between p° and /o(980) or a broad scalar cr(400 — 1200) are estimated to be small. The 
orbital angular momentum for p'^Tr" (p^vr"'") is one, while for /o(980)7r'^ or air^ (/o(980)7r'^ or air'^) 
it is zero. Therefore the two wave functions are orthogonal. The interference term vanishes when 
integrated over the whole space. As a cross check, MC samples with interference effects are made 
from non-resonant S+ — tt'^'tt^tt"'" and B^ — 7r"^7r~7r° Monte-Carlo using a reweighting technique. 
The full selection is then applied. The relative phase is chosen to maximize interference. Small 
effects are observed, as expected. 

Table 5 summarizes the various sources contributing to the systematic errors in the branching 
fractions. The dominant systematic errors are due to the fit procedure (imperfection in likelihood 
model) and the uncertainties in the B background model. Table 6 summarizes the possible sources 
contributing to the systematic errors in the charge asymmetries. 
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Table 5: Results and breakdown of systematic errors for the branching ratios measurements. 



Signal yields and efficiencies 





B+ ^ p+ttO 


B+ ^ p''tt+ 


5° ^ p^ttO 


Corrected signal yield 


170.8 


232.5 


15.6 


Corrected esignal 


17.5 ± 0.1% 


28.3 ± 0.1% 


20.0 ±0.1% 


Statistical error on signal yield 


OS 7 
Zo. 1 


26.4 


11.7 




Yield 


systematics (absolute) 


Tfiib 0.016 ps 


n/a 


n/a 


0.1 


At resolution model 


n/ a 


n/a 


0.3 


B tagging 


3.7 


n/a 


0.9 


Fraction of misreconstructed signal 


3.0 


1.3 


0.1 


AE PDF 


8.5 


0.7 


1.0 


ruES PDF 


2.5 


1.6 


1.5 


NN PDF 


3.0 


3.2 


2.6 


B backgrounds 


11.2 


2.3 


+3.2 
-3.9 


Fitting procedure 


14.4 


8.2 


6.3 




21.0 


9.3 


8.1 




Relative efficiency and scaling 


; systematics 


Tracking efficiency correction 


0.8% 


2.4% 


1.6% 


PID for tracks 


1.7% 


5.2% 


4.0% 


Neutral correction 


10.2% 


n/a 


5.1% 


AE cut efficiency 


2.6% 


1.0% 


0.1% 


rUES cut efficiency 


0.0% 


0.0% 


0.0% 


p cut efficiency 


1.5% 


0.0% 


0.0% 


NN cut efficiency 


4.0% 


4.0% 


1.0% 


DIRC photons cut for bachelor 7r+ 


n/a 


1.0% 


n/a 


B{p TT+TT") 


1.6% 


1.6% 


1.6% 


N{BB) 


1.1% 


1.1% 


1.1% 


Sub-total 


11.6% 


7.2% 


7.0% 


Total systematic error 


16.9% 


8.3% 


52.5% 


Branching ratio [xlO~^] 


11.0 ± 1.9 ± 1.9 


9.3 ±1.0 ±0.8 


0.9 ±0.7 ±0.5 



20 



Table 6: Results and breakdown of systematic errors for the measurements of charge asymmetries 
in B+ p+tt'^ and B+ p%+ . 





B+ ^ /9+7r0 


B+ ^ p07r+ 


Acp 


0.23 


-0.17 


Statistical error on Acp 


0.16 


0.11 




Fit systematics 


B tagging 


0.004 


n/a 


Fraction of misreconstructed signal 


0.003 


0.0012 


AE PDF 


0.034 


0.0004 


THeS rUr 


0.003 


0.0016 


NN PDF 


0.00 


0.0028 


B backgrounds 


0.050 


0.022 


Detector charge asymmetry 


0.01 


0.009 


Total systematic error 


0.06 


0.02 


Acp 


0.23 ± 0.16 ± 0.06 


-0.17 ±0.11 ±0.02 



9 Summary 

We have presented preliminary measurements of branching fractions and CP- violating charge asym- 
metries in S+ — > /3'*"7r° and 5+ — > p^Tr"*" decays, and a search for the color-suppressed decay 
^0 _^ pOyrO. The data sample used in the analyses consists of 89 x 10^ BB pairs. We find a 
branching fraction for B~^ that is consistent with previous measurements [12, 13]. We 

observe a signal for B~^ p'^Tr^ with a significance of 6.6(7, and set an upper limit for B^ — > p^ir^. 
We do not observe evidence for direct CP violation. 
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